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ABSTRACT 



OO ' We analyze three sets of infrared star counts in the inner ~ 0.5 pc of the Galactic 

Center. We perform statistical tests on the star counts and model in detail the 
extinction field and the effects of dwarf-giant collisions on the luminosity function. We 
find that both the star counts and the depletion of the brightest stars in the inner 
'sj" . ~ 0.05 pc can be explained by a n oc r~° stellar cusp with a in the range 3/2 to 7/4, 

^ ! in which the envelopes of the brightest giants are destroyed by stellar collisions. Such 

O ! a cusp is consistent with the Bahcall-Wolf solution for the distribution of stars that 

have undergone two-body relaxation around a black hole. We show that systematic 
uncertainties due to variable extinction and unrelaxed stars are probably small, but 
CL|' deeper star counts are required to confirm these results. We estimate that the tidal 

Q ■ disruption rate of cusp stars by the black hole is few x 10^^ yr~^. 

H ■ 

^ I Subject headings: Galaxy: center — Galaxy: kinematics and dynamics — 

' Galaxy: stellar content — infrared: stars 



1. Introduction 

The Galactic Center (GC) offers a unique opportunity to study the dynamical effects of a 
super-massive black hole (BH) on the stellar population in its immediate vicinity. Dynamical 
models of the evolution of such a system generically predict the formation of a stellar cusp 



with a characteristic slope that reflects its formation history (Bahcall & Wolf 1976, 1977; Young 



1980; Lee & Goodman 1989| ; Quinlan, Hernquist & Sigurdsson |1995| ). Although the existence 



of super-massive BHs in the center of normal galaxies now appears to be the rule rather than 



the exception (Magorrian et al. |1998D , the observed central light distributions do not seem to 



be correlated with the mass of the central black hole (Kormendy & Richstone 1995| ; but see 



Magorrian et al. 1998). It appears that even on the smallest resolved scales of about 0.1", the 
light distribution still reflects the overall galactic dynamics, and that the effects of the BH on the 
spatial distribution of the stars may be confined to yet smaller scales. Even with better resolution, 
it will be very hard to observe individual stars in the dense galactic cores, and the properties 
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of the stellar population will have to be inferred from the projected integrated light, which is 
completely dominated by giants and supergiants. These rare stars may not be representative of the 
general stellar population. Furthermore, their light-emitting envelopes could be affected by various 
processes in the extreme environment inside a BH cusp. The situation is further complicated 
by the fact that galactic centers often display cusps in the light distribution on scales that are 
much too large to be related to the presence of a super-massive BH in the core. Obviously, other 
dynamical processes can also give rise to stellar cusps, so it is not possible to claim that a central 
stellar cusp must be related to the presence of aBH (Phinney |l989| ). 

This is to be contrasted with the situation in the GC, where precise measurements of the 
apparent motions and radial velocities of individual stars have already made it possible to weigh 



the central BH dynamically, independently of the light distribution (Genzel et al. |1997| ; Ghez et 
al. 1998|). This makes it legitimate to interpret the stellar distribution very near the center in the 



context of BH formation scenarios. These observations also yielded deep infrared star counts, which 
extend the stellar census to stars far less luminous and much more numerous than supergiants. 
However, the step from the observed number counts to the actual stellar distribution is far from 
trivial. Observations show that the stellar population of GC is a mixture of old and young stars, 
which are probably the products of episodic star-formation throughout the GC history (Genzel, 
Hollenbach & Townes |1994 ; Serabyn & Morris |1996D . Some of the young and very luminous stars 
are not dynamically relaxed (Krabbe et al. |1995| ; Genzel et al. 1996 ). Spatial variations in the 
spectral properties of the stellar population (Sellgren et al. |1990| ; Krabbe et al. 1995 ; Genzel et 
al. 1996|) , and possible spatial gradients in the extinction (Blum et al. 1996 ; Davidge 1998| ), could 
further complicate the interpretation of the star counts. The necessary first step towards testing 
dynamical models in the inner GC is to establish what is the stellar distribution there. 

The evidence for a stellar cusp within the cusp radius of the GC (roughly, the radius enclosing 
a stellar mass comparable to the BH mass, which is ~ 3 pc in the GC) has so far been inconclusive 
or negative. The mass models that are currently adopted in stellar proper motion studies of the 
GC do not have an inner stellar cusp, but assume a smoothed isothermal distribution with a core 
radius of 0.4 pc (1" = 0.04 pc in the GC for a solar galactocentric distance of 8 kpc (Carney et al. 



1995 ). In this work we re-examine the question of whether there is a stellar cusp in the inner 10" 
of the GC. 

This study is motivated by several recent developments. Progress in extending infrared 



observations of the GC to fainter magnitudes and higher angular resolutions (Genzel et al. |1997| ; 
Ghez et al. 1998| ) and to larger samples with multi-color information (Blum et al. 1996; Davidge 



et al. 1997 ; Davidge 1998 ) yielded a wealth of information on the distribution of stars in the GC. 
In addition, these deep infrared counts in the GC, together with infrared and optical star counts 
in Baade's window (Tiede et al. 1995| ; Holtzman et al. 1998| ), now make it possible to reliably 
model the infrared stellar luminosity function in the GC over a large luminosity range (Alexander 
& Sternberg 1999| ). However, the high angular resolution studies are primarily focused on tracking 
stellar proper motions for weighing the BH, while the multi-color studies are focused on studying 
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the spectral properties of the stellar population. Consequently, the star counts, which are a 
by-product of these studies, have not yet been analyzed in detail, and the additional information 
supplied by the luminosity function has not been taken into account. Finally, recent high resolution 
deep imaging of the central few arcseconds revealed an apparent over-density of faint stars in the 
inner l" relative to the immediate surroundings, the so-called "Sgr A* cluster" (Genzel et al. 1997] ; 
Ghez et al. 199^ ). This raises the question whether a distinct stellar cluster with an atypically 
faint stellar population could have possibly formed in the inner 1" or else migrated there from 
farther away. 



It this study we analyze the observed stellar distribution in three published data sets (§p.3| 



Keck data from Ghez et al. (|1998D , OSIRIS data from Blum et al. (|1996D and SHARP+3D data 
from Genzel et al. ( |1996|) and Eckart & Genzel ( |1997l ). As we will argue below, these data point 
to a simpler explanation for the Sgr A* cluster, namely that it is the tip of an underlying stellar 
cusp that smoothly rises throughout the inner 10", and that the absence of luminous stars there 
is due to a combination of envelope-destroying stellar collisions and projection effects rather than 
an intrinsic change in the properties of the stellar population. We will further show that the 
shape of the cusp agrees with that predicted by the Bahcall-Wolf solution for the distribution of 
intermediate mass stars in a multi-mass system that has undergone two-body relaxation around a 



BH (Bahcall k Wolf [1976| , |1977| ; Murphy, Cohn & Durisen |1991D , as appears to be the case in the 
inner GC. 

This paper is organized as follows. In we describe the relevant properties of stars in the 
GC and present the data that are used in our analysis. In ^ we describe the stellar density 
distributions that we use to model the data and present the results of the statistical analysis of the 
counts. Some anomalies in the counts are further discussed in Appendix |^. In §^ we investigate 
the effects of extinction on the star counts. Readers who are mainly interested in the dynamical 
results of this work may prefer to skip this section and proceed to where we discuss the effects 
of stellar collisions on the depletion of the brightest stars. The orbit ally- averaged collision rate is 
derived in Appendix |B|. We discuss the results in ^ and summarize our conclusions in 



2. Stars in the Galactic Center 

2.1. The stellar population and stellar distribution 

The stellar population in the GC is a mixture of old Galactic Bulge stars and younger stars 
that were produced in various star-formation episodes during the lifetime of the Galaxy (Genzel, 
Hollenbach Sz Townes 1994 ; Serabyn & Morris 1996|) . About 80% of the stars observed in the 



central parsec down to i^T = 14.5™ are identified as late-type giants and supergiants by their 
pronounced CO absorption. Late-type stars at the magnitude range 11™ ^ K ^ 14.5™ are K 
and M giants, which correspond to a mass range of l-TM© and ages > few x 10^ yr. Late-type 
stars at the magnitude range 9™ ^ K ^ 12™ are probably asymptotic giant branch (AGB) stars 
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in the mass range 2-8Mq and ages of 1-few x 10^ yr. The most luminous late-type star, IRS7 
{K = 6.5™), is a massive 15-25Mq supergiant with an age of no more than 10'' yr. The late-type 
giants are dynamically relaxed and follow the galactic rotation (Genzel at al 1996 ). 

The remaining ~ 20% of the K < 14.5™ stars in the central parsec appear to be recently 
formed stars, either exhibiting strong Hel/Hl emission or a steep featureless reddened continuum 
(Krabbe et al. 1995; Genzel et al. 1996| ). Two extended red objects, which may be recently 
formed stars still embedded in a dust shell, are also observed (Ott, Eckart & Genzel 1998|) . The 
Hel/Hl stars are massive, hot and very luminous stars with huge mass loss rates, which are 
probably undergoing a short lived phase lasting less than 10^ yr, on their way to becoming classical 
Wolf-Rayet stars. In this case, their progenitors must be ^ 20Mq stars with main-sequence 
lifetimes shorter than 10'^ yr. The ~ 20 Hel/Hl stars in the inner few arcseconds supply more than 
50% of the total 2 ^m flux there. The fainter stars with the featureless reddened continua could 
be main-sequence O stars. The four stars in the Sgr A* cluster (inner 1") for which spectra were 
obtained appear to be such early-type stars (Genzel et al. 1997| ). Krabbe et al. (1995) modeled 
the Hel/Hl stars and the late-type supergiants in the central parsec by a SOOOMq starburst 
that occurred between 3 and 8 x 10^ yr ago. This is a very small fraction of the total stellar 
mass that is enclosed in the central parsec, ~ 1.5 x W^Mq. These recently formed stars are not 
dynamically relaxed, as is demonstrated by the radial velocity field of the Hel/Hl stars in the inner 
~ 10", which clearly indicates that the stars are counter-rotating at a constant rotation velocity 
Viot = 120kms~^ with respect to the galactic rotation of Vrot = 24kms~^. There are additional 
indications that the stellar population varies with position and is not quite well-mixed. The depth 
of the CO absorption feature in the integrated light decreases in the inner ~ 10" (Sellgren et 
al. 1990| ). This is probably due to the joint effects of a depletion of the luminous {K < 10.5™) 
late-type stars in the inner 5" and the dilution of the absorption feature by the concentration of 
luminous early- type stars there (Krabbe et al. 1995| ). 

The basic premise of the attempts to measure the global properties of the stellar distribution 
in the inner GC is that there is an underlying dynamically relaxed distribution, which includes 
most of the mass and most of the stars, and which is traced by the old faint stars. In contrast, 
the young stars and bright giants, which dominate the light, are less reliable tracers because they 
could be unrelaxed, or subject to statistical fluctuations due to their small numbers, or susceptible 
to environmental effects such as stellar collisions or photoionization (Sellgren et al. 1990| ). 

Our current understanding of the stellar distribution in the GC is incomplete. The light 
distribution approximately follows a power-law from as far away as ~ 100 pc, where it cannot 
be directly related to the central super-massive BH, to ~ 1 pc from Sgr A* (Serabyn & Morris 



1996). Previous attempts to characterize the distribution in the inner parsec concentrated on 



identifying the core radius, rc , where the surface density falls to half its central value. The value 
of Tc is much debated in the literature, and ranges from 0.05 pc (Allen 1983| ) to ~ 1 pc (Rieke & 
Rieke 198^ ). It appears that the differences between the various estimates depend on whether the 
light distribution or the number counts are used and whether or not the bright stars are included 
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in the fits. 

The most recent estimates for the mass distribution in the inner GC are dynamical and are 
based on the stellar apparent motions and radial motions. Genzel et al. ( 1997 ) and Ghez et al. 
(1998) fit the enclosed mass estimates to a flattened isothermal distribution of the form 



P^'^ = 1 + 3(r/r.)2 ' 

where r is the 3-dimensional distance from the BH, pQ = 4 x IO^Mqpc"^ and Tc = 0.4 pc. Because 
the BH mass completely dominates the velocity field at r <C r-c, and the stellar mass begins to 
affect the stellar velocities only at r ~ Tc, this mass distribution measures only the integrated 
stellar mass up to r ~ Vc, and is insensitive to the distribution of the stars within the core volume. 
Therefore, the good fit of the velocity field to this mass model does not necessarily preclude the 
existence of a stellar cusp in the core. 



2.2. The K luminosity function and mass function 



The construction of the K luminosity function (KLF) and mass function models for the GC 
is discussed in detail by Alexander & Sternberg ( 1999| ). Briefly, the observed KLF in the GC 
(Blum et al. |1996| : Davidge et al. |1997D can be smoothly joined to that in the Galactic Bulge, as 



observed through Baade's window (Tiede et al. |1995D , and further extended to lower luminosities 
by transforming deep optical counts in the bulge (Holtzman et al. |1998D to the i^-band. The 
resulting KLF (fraction of stars per K-magnitude) has the form df /dK oc lO''^ with b = 0.35, and 
extends from brightest observed stars in the GC at Ki = 8™ down to the main-sequence turnoff 
at K2 ~ 21.5™. Stars below the main-sequence turnoff at ~ I Mq have negligible contribution to 
the KLF. By normalizing the KLF to the observed star counts in the GC, it can be verified that it 
reproduces the observed mass-to- light ratio in the GC. This empirical KLF and mass-to-light ratio 
can also be reproduced remarkably well by a population synthesis model that assumes a continuous 
star-formation history with a Miller- Scalo (Miller & Scalo 1979; Scalo 1986| ) initial mass function 
(IMF) with a mass range of O.IMq to 120Mq. The IMF can be then used to follow the evolution 
of the stellar mass function over time (§^). 

We adopt these luminosity and mass functions as a volume averaged description of the 
stellar population in the inner parsecs of the GC and assume that they adequately describe the 
"well-mixed" population. 



2.3. The observed star counts 
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2.3.1. The Keck data set 

Ghez et al. (|199^ ) obtained diffraction limited (0.05" resolution) K-band images of the 
inner 5" x 5" with the Keck interferometer, down to a detection threshold of Kq = 17™. Using 
conservative selection criteria, they extracted from the field 90 stars for which the angular motions 
could be reliably measured and tabulated their magnitude and projected positions and velocities 
(Figure |l]). As the field is slightly off-center with respect to Sgr A*, only the inner projected ~ 2.5" 
are covered over the full 2tt radians. The original criterion for inclusion in the data set was the 
reliability of measuring angular motion, rather than completeness in magnitude. Consequently, 
the counts are complete only down to ~ 14™-15™, beyond where the completeness falls off very 
rapidly (A. Ghez, private comm.) to a level of ^ 1% in the range ~ 16™-17™. The incompleteness 
is not uniform over the field, as it depends on the distribution of the bright stars. In our analysis 
we consider only the 64 stars in the inner 2.5", of which 10 are identified as early type stars. We 
assume here that the remaining 54 spectrally unidentified stars are late type stars. 



2.3.2. The OSIRIS data set 

Blum et al. ( |1996D imaged the inner 2' of the GC with the Ohio State Infrared Imager and 
Spectrometer (OSIRIS) in the J, H and K bands (~ 1" resolution) for the purpose of measuring 
the luminosity function and color-magnitude diagram of the stellar population. They detected 
a total of ~ 2000 stars and determined by simulations that the entire sample is complete down 
to -fC ~ 12™. Of this large sample, they tabulated the K magnitude, projected positions and K 
extinction coefficient Ak for a sub-sample of 147 stars, which span a range of K = 6.4™-12.4™. 
The stars included in the sub-sample are those brighter than K = 10.5™, those for which Blum 
et al. obtained 2^m spectra, or those with IRS numbers. Thus, this sub-sample is guaranteed 
to be complete only down to K = 10.5™. The distribution of the stars in the field is very sparse 
beyond projected radius p ~ 10"-15" (Figure §). We limit our analysis to the 50 late type stars in 
the inner 13" in order not to be affected too much by the Poisson fluctuation in the counts or by 
molecular gas and dust that may extend inward from the circumnuclear ring (Genzel, Hollenbach 



Townes 1994). 



2.3.3. The 3D and SHARP data set 

Genzel et al. ( tl996| ) used the 3D imaging spectrometer to identify spectrally 223 stars in 
the central parsec and separate them into early and late-type stars. The stars were assigned to 
approximate magnitude classes: bright {K < 10.5™), intermediate (10.5™ < K < 12™) and faint 
(12™ < K). Because the 3D counts are incomplete in the inner 5" due to source confusion, they 
augmented the counts with high resolution SHARP camera observations of < 13.5™ stars in the 
inner 10" (Eckart & Genzel p^97| ), but did not include the 11 very faint {K > 14™) Sgr A* cluster 
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stars (S'l-S'll) in the inner ~ 0.5" (A. Eckart, private comm.). With the faintest stars excluded 
from the sample, the composite data set is of roughly uniform depth over the central ~ 20". 

Genzel et al. ( |1996| ) binned the composite data set and fitted a flattened r~^-^ power-law 
distribution (analogous to Eq. ||) to the surface number density and obtained a core radius of 
rc = 0.29^0,11 pc (rc = 7.4" llls)^- We reconstruct here the inner ~ 20" of the the unbinned 
composite data set for from the tabulated 3D and SHARP stellar positions in order to check 
whether other density models can fit these star counts as well. Our composite data set includes in 
total 213 stars, of which 198 are late-type 3D stars from Genzel et al. ( 1996 ) and 15 are spectrally 
unidentified SHARP stars from Eckart & Genzel (1997). The surface distribution of the 213 stars 
is shown in Fig. ^. To avoid edge-of-field effects, we consider in our analysis below only the 107 
late-type stars in the inner 13", which are composed of 92 3D late-type stars and the 15 SHARP 
stars. 



3. Statistical analysis of the star counts 

3.1. Stellar density distribution models 

The main purpose of this work is to quantify the empirical stellar density distribution in the 
GC, but theoretical expectations guide us in the choice of models for describing the distribution. 
Generally, the effect of a super-massive BH on the stellar population within the cusp radius will 
depend on whether the relaxation time due to two-body stellar encounters is longer or shorter than 
the system's age. When two-body relaxation is unimportant, the black hole may grow adiabatically, 
and the resulting stellar distribution will depend on the initial conditions. A spherically symmetric, 
non-rotating system that is initially isothermal will develop a r~^/^ cusp (Young [198Q|) . Other 
initial conditions will result in a variety of density distributions, in some cases as steep as r~^/^ 
(Lee & Goodman 1989| ; Quinlan, Hernquist &: Sigurdsson 1995| ). 



When two-body encounters are effective, the final configuration will not depend on the initial 
one. Very close to the BH, where inelastic stellar collisions dominate, a r~^/^ cusp of marginally 



bound stars will form (Murphy, Cohn & Durisen 1991). At larger radii, where collisions are elastic 



a single mass population will settle into a r"''/'' cusp (Bahcall &: Wolf 1976). In a realistic stellar 
population with a spectrum of masses, mass segregation will take place, and the massive stars 
will be more concentrated towards the center than the low-mass stars. Bahcall Sz Wolf ( 1977|) 



investigated the two-body relaxation of a two-component stellar population around a massive BH. 
They found that the two components settle into different power-law distribution with power-law 



^ This is a somewhat smaller value than rc — 10" that was adopted in their subsequent works as a compromise 
between the different core radii suggested in the literature (R. Genzel, private comm.). 
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indices given by 



mi 3 



Murphy, Cohn & Durisen (|l99lD showed that Eq. g holds also for the spatial distribution of a 
multi-mass population, when m2 is taken to be the maximal mass cutoff, so that 3/2 < ai < 7/4. 
It is straightforward to verify that over a limited range in r, the total density distribution that will 
be observed will still be an approximate power-law with an index in this range. 

The local two-body relaxation time at the core radius of the dynamical mass model (Eq. |l]) 
can be estimated by (Spitzer & Hart 1971J ) 



' G2(m)plnA 

where A ~ M,/ (m) and we assume (m) = 0.6 Mq, a BH mass of M, = 2.6 X 10*^ Mq (Genzel et al 



19971 ; Ghez et al. |1998|) and p = 10^ Mopc'^, a = 117kms-^at rc = 0.4 pc (Genzel et al. |1996| ). 



Since the relaxation time is shorter than the age of the Galaxy, the theoretical expectation is that 
the observed giants, which have an intermediate mass, will have a power-law distribution with 
an index in the range 3/2 < a < 7/4. However, we will avoid here any theoretical prejudice and 
investigate power-law distributions with indices that span the full range of the various theoretical 
scenarios from a = to 5/2. 

We consider two spherically symmetric stellar density models. The first is the flattened 
isothermal model (Eq. ||), where we assume a constant ratio between the mass density and the 
stellar number density 



n(r) 



no 



l + 3(r/rc)2' 



(4) 



and where no is the central number density and Tc 
surface density is 



10" (Genzel et al. 1996). The corresponding 



where p is the projected distance from the center. The second class of models are broken 
power-laws, 



ro < r < 
r>n, (3>1 



(6) 
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where Ub is the density at the break at ri, , and ro is the closest distance to the BH where stars can 
exist. Although rf, is not the core radius, we will loosely refer to the region interior to as the 
core. The corresponding surface density is 



x{H{a,^)(^f''-H(p,^)(f-y'] ro<p<rb (7) 



2' 



r^B i^)p^~^ p>rb 



where B is the beta function and H is defined by the hypergeometric function 



Hix,y) =2Fi(i 



3.2. The cumulative distribution function 

Figure § shows the cumulative distribution functions (DFs) for the late-type stars in the 
three data sets (more precisely, all the observed stars that were not positively identified as Hei/Hi 
stars, early- type stars or dust embedded stars). The cumulative DFs indicate that neither a 
flattened isothermal model (Eq. ^) nor a singular isothermal model n oc provide a good 
description of the observed star counts in the central ~ 10". This impression is supported by 
the Kolmogorov-Smirnov (K-S) test results (presented as acceptance probabilities in Table [I]), 
although the K-S scores do not differentiate strongly between the different models. Generally, the 
K-S scores favor cusps with slopes in the range a ~ 3/2 to 7/4. The discrepancy between the 
OSIRIS data and the flattened isothermal model is the smallest of all data sets, although this is 
not the best-fitting model. This probably reflects the fact that the OSIRIS data set is less deep, 
and is therefore more sensitive to the decrease in the number of luminous late-type giants in the 
central few arcseconds. 

The statistical significance of the concentration of stars at p ~ 2", which is seen in the 
cumulative distribution of the Keck data set, is discussed in appendix 



3.3. Maximum Likelihood analysis 

The trends that are seen in the cumulative distributions of the star counts (Figure I) and 
the theoretical expectations (§^?l|) suggest fitting power-law density distributions to the star 
counts. In order to avoid the arbitrariness of binning, we perform a Maximum Likelihood (ML) 
analysis to find the best fitting power-law index a, assuming a broken power-law distribution 
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with fixed values (3 = 2 and ri, = 10" (we verified that the results are not very sensitive to the 
exact shape of the outer distribution). Figure || shows the likelihood curves that were obtained 
by using Eq. ^ to calculate the relative probabilities of the observed stellar positions^. As can 
be expected, the ML best-fit parameters are not exactly the same as those derived from the 
Kolmogorov-Smirnov analysis, since the two methods are sensitive to different ranges of the data 
(the K-S to the mid-range of the projected distance and the ML to the inner, high-probability 
region). Nevertheless, the results are qualitatively similar. The likelihood curves for all the samples 
peak in the range a ~ 1.5 to 1.75, irrespective of whether the they are based on stars in the inner 
13", or only the inner 2.5". It is reassuring that this range is consistent with the ML estimate for 
the entire Keck data set, as this suggests that the contamination by the ~ 20% of the young stars 
does not introduce a large bias. A flat core (a = 0) is ruled out only at at the 2a level by the Keck 
and OSIRIS data sets, but is completely incompatible with the 3D+SHARP data. 

We also performed a ML analysis to find the best fitting core radius Tc, assuming a flattened 
isothermal distribution. The results are qualitatively similar to those of the broken power-law. In 
all data sets the likelihood curves peak at rc < 2.5", and while a 10" core radius is still compatible 
with the Keck and OSIRIS data sets at the ~ 1.5cr level, it is ruled out by the 3D+SHARP data 
at the ~ 2.5(T level. 

The best fit core radius we derive from our composite 3D+SHARP data set for the late-type 
stars (rc = 1.5") is much smaller than the one derived by Genzel et al. ( |1996| ) (rc = 7.4"l^;g). 
The reasons for this disparity are unclear. The most significant difference between the two analysis 
methods appears to be the issue of binning. The binned star counts that were used in the fit 



of Genzel et al. (|1996D could lead to the suppression of a steep cusp in sparse counts. Our initial 
attempts to fit density models to binned 3D+SHARP data indicated that the best-fit density in 
the innermost regions fluctuates considerably with small changes in the bin sizes or bin centers. 
A more detailed comparison of the unbinned ML and the binned fitting procedures, which is 
outside the scope of this work, will be needed to resolve this issue. 



4. Extinction in the inner GC 

A gradient in the extinction field could bias the star-counts because fewer stars would be 
observed where the extinction is higher. For example, a central concentration of dust could make a 
stellar cusp appear as a flat core, and conversely, a dust ring could make a flat core appear as a 
cusp. In this section we use observations of the extinction in the GC to estimate the magnitude of 
such a bias. 



^ For each density model, we normalize the probability to the projected distance interval between the innermost 
and outermost observed stars in the sample. This is necessary for removing biases due to the finite size of the field, 
and biases due to divergent behavior of the density at p = 0, where the measured positions are not much larger than 
the angular resolution and the ~ O.l" uncertainty in the position of the dynamic center (Ghez et al 1998). 
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Davidge ( 1998 ) used J, H, and K observations of stars in the inner 3' to estimate the radial 
dependence of the K-band extinction coefficient, Ax, averaged over 20"-wide bins. Davidge found 
that {Ak) systematically increases from 2.8™ at ~ 100" to 3.2™ at ~ 30", and then decreases again 



at the inner 20" to 2.9™ (AAi^- ~ —0.02™ per arcsec). Blum et al. ( |1996 ) made a similar study of 



the inner 2', and obtained an average central extinction value of (Ak) ~ 3.3™. Eckart et al. ( 1995 ) 
found no significant color gradient over the Sgr A* cluster in the inner l" . Figure ^ shows the 
distribution of stars in the Ak — p plane in the inner 25" of the GC (Blum et al. tL996|) . A small 



fraction of the stars, all concentrated in the inner 10", are heavily reddened. One of those (IRS 21) 
was identified by Ott, Eckart & Genzel ( 1998|) as an extended and highly polarized object, which 



may be a recently formed star still embedded in a dust shell. A second such star (IRS IW) was 
also identified in the inner 10", but is not part of the OSIRIS data set. 

The observed total surface density, Sq, and the observed differential surface density, dTio/dAx, 
reflect the 3-dimensional distribution of the dust. The ii'-extinction suffered by a star at projected 
distance p and at point z along the line-of-sight can be decomposed into two parts. The first 
term, Ai{p), corresponds to extinction by a foreground dust screen, and the second term, A2{p,z), 
corresponds to that by dust inside the GC (Figure]^. The total extinction, Ak, is then 



PCX) 

Ak{p, z) = Ai{p) + A2{p, z) =Ai{p) + D / d{p, z')dz' , 

J z 



(9) 



where d is the dust density and D is a factor that translates the dust column to extinction. 
Assuming that the unreddened K luminosity function (KLF) is constant over the volume of interest 
and assuming cylindrical symmetry around the line of sight, the distribution of stars in the p — Ak 
plane is given by 2TTp{dT,Q/dAK), where 



and Ko is the detection limit. The observed total surface density So will be different from the true 
one since 



/r-Ko-AK(p,z) Jf 
n{p,z)dz -^dK. (11) 

It is important to emphasize that the measured quantity {Ak{p))-, which is averaged over the 
observed stars, is not directly related to the average expressed by Eq. |ll|, since regions of high 
extinction will be under-represented in the counts, and since intrinsically very luminous stars, 
which can be observed even through high extinction, are very rare. Therefore, the relatively small 
measured value {Ak) ~ 3.5 does not necessarily preclude the possibility that the innermost GC 
is hidden by strong extinction. In order to test this, it is necessary to consider the 3-dimensional 
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distribution of the dust. Because it is unclear what completeness corrections have to be applied to 
the data presented in Figure ^, we will restrict ourselves here to qualitative comparisons of some 
simple dust geometries with the observed distribution of stars in the p — Ax plane. 

We consider first the case where the dust is concentrated in a foreground screen, so that 
Ak{p,z) = Ai{p). If Ai is a constant, or varies randomly over small enough angular scales, Sq 
will simply be proportional to S, at least on average, and the extinction will not bias the observed 
density distribution. This requires us to interpret the few observed highly reddened stars as 
intrinsically reddened objects that do not trace an extended dust distribution. A linear regression 
of Ax on p (Figure does not show any strong systematic trends in the radial behavior of Ax 
in the inner 20" and is thus consistent with a foreground dust screen, where the scatter is due to 
random fluctuations in the extinction, or measurement scatter, or both. The differential extinction 
is /S.Ax = +0.01™ per arcsec based on the the inner 10", and IS.Ax = —0.03™ per arcsec based on 
the inner 20", and in both cases /^Ax = is within the la errors. 

The bias due to a non-zero trend of this order will not be large. If a power-law density 
distribution is estimated from the counts at points p\ and P2, then the difference between the 
measured power-law index and the true one is 



^^^ Km-n)^Ax 

loglo(P2/Pl) 

which is only Aa < 0.1 for /S.Ax = —0.03™ per arcsec and pi = 1", p2 = 10". 

Next, we consider the general case, where in addition to a foreground dust screen the BH and 
the stars near it are embedded in a dust cloud. In particular, we wish to address the possibility 
that the highly reddened stars in the central 10" are those few stars that are luminous enough 
to be observed even from inside the dust cloud. We assume a simple configuration where Ai is 
constant and the dust is distributed in a spherical cloud with a density profile 

^ = 1 + ^\r/RY ' (1') 

where R is the core radius of the dust cloud (Figure |^ . Figure |8| shows the relative probability 
density for observing stars in the p — Ax plane, as calculated from Eqs. |9| and |l^, for a n oc 
stellar distribution and a flattened isothermal distribution with a core radius of 10" (Eq. ^). These 
two stellar distributions were chosen to represent the two cases of a sharp density cusp and a flat 
core. We set the detection threshold at Kq = 13™, corresponding to the OSIRIS data threshold, 
and find by trial and error that the parameter values D = 0.15, R = 30" and Ai = 2.5™, which 
correspond to a maximal extinction of 10.7™, produce a large enough spread in Ax at ~ 10". Dust 
models with lower values of R and D fail to reproduce the highly reddened stars because they do 
not have a high enough dust column, while in models with higher values of R and D the extinction 
is so effective that the central 10" is almost completely obscured. Figure ^ show that the models 
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differ from the observations in several significant aspects. Both predict a strong concentration of 
stars near and above Ai, which increases with p as the extinction decreases. Such an increase is 
not observed in the data. Both also predict a nearly constant distribution of high-A^ stars from 
p ~ 5" outwards (due to the fact that both the dust and stellar columns fall together), unlike the 
observed distribution, which is confined to the inner 10". The flat core model predicts very few 
high- Ax stars in the inner ~ 5", since the flat stellar distribution does not provide enough high 
luminosity stars that can be observed through the dust cloud. This is to be contrasted with the 
dusty stellar cusp model, which does predict a sparsely populated region in the p — Ak plane at 
p ^ 10" and Ak ^ 7™, approximately as is observed. However, the stellar cusp model also predicts 
that a strong concentration of highly reddened stars at p ^ 2", contrary to the observations. We 
conclude that distributions that are dusty enough to account for the highly reddened stars inside 
10" are not favored by the observed distribution of stars in the p — Ak plane. 

However, the p — Ak distribution in itself cannot rule out a combination of a moderately dusty 
cloud together with a few intrinsically reddened stars, and furthermore, our conclusions depend on 
the assumed stellar DF and the KLF. The few spectra that were obtained for stars deep in the 
potential well can provide additional information that is independent of these assumptions. Such 



stars can be identified by their high velocities. Genzel et al. (|1996|) obtained spectra of 4 stars in 
the Sgr A* cluster (SI at p = 0.11" with projected velocity v± = 1556 km s~^, S2 at p = 0.17" with 
v± = 691 km s-\ S8 at p = 0.32" with v± = 1086 km s'^, and Sll at p = 0.54" with v± = 960 



km s ). Assuming that the stars are bound to the M, = 2.6 x 10^ Mq BH (Genzel et al. |1997| : 



Ghez et al. |1998D , these velocities translate into maximal 3-dimensional distances from the center 
of rmax = 2GM,/vl = 0.24", 1.2", 0.48", and 0.62", respectively. The four spectra are featureless 
and appear to be consistent with those of moderately reddened {Ak ~ 3™) early-type stars in the 



central cluster (Genzel et al. |1996D . This suggests that there is no significant concentration of dust 
around the BH. 

We conclude that spherical dust distributions that are dusty enough to account for the 
luminous, highly reddened stars in the inner 10" require a stellar cusp to provide enough of such 
stars and are therefore inconsistent with the absence of a strong concentration of highly reddened 
stars in the inner few arcseconds. In addition, the absence of a color gradient in the inner 1" and 
observations of moderately reddened stars deep in the potential well further argue against a central 
dust cloud. On the other hand, the observed p — Ak distribution is consistent with extinction 
by a uniform foreground dust screen, provided that some intrinsically reddened stars exist in the 
inner 10". Two such stars have in fact been detected there. A foreground dust screen will not 
introduce a large systematic bias in the determination of the effective power-law index of the stellar 
distribution. 
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5. Stellar collisions and the depletion of luminous stars 



The problem of the incompleteness of the counts could be circumvented, in principle, by 
restricting the analysis to the brightest late-type stars, which should directly reflect the projected 
stellar density. However, stars brighter than K = 13.5™ are not observed in the inner 1" of the GC 
in the Keck data set, although the counts are complete at least down to K = 14™. Figure ^ shows 
the distribution of stellar magnitude against projected radius. Apart from the counter-rotating 
Hei/Hi stars, which clearly stand out in the p — K plane, there is a gradual trend towards a fainter 
maximal luminosity in the inner 2". The possibility that this is due to a gradient in the extinction 
is not supported by the analysis in §^ 

Another interesting possibility is that the luminous stars have been destroyed by stellar 
collisions. Stellar collisions were proposed as the mechanism for destroying the outer envelopes 
of late-type giants^ (Lacy, Townes &: Hollenbach 1982 ; Phinney 1989 ) as an explanation for the 
depletion of the CO absorption in the integrated light (Sellgren 1990| ) or for the depletion of the 



bright late-giants in the counts in the inner 5" of the GC (Genzel et al. 1996 ). However, the 
dynamical timescale in the envelope of a giant star is of the order of a year and the thermal 
timescale is of the order of several decades, so that a collision in which a star "punches a hole" in 
the giant's envelope will probably not have a lasting effect. Studies of collisions between dwarf 
stars and giants indicate that a significant fraction of the collisions can result in a permanent 
destruction of the envelope if the ratio between the impact parameter and the giant star radius, Xe, 
is small enough, ^ 0.5. The envelope is destroyed either by the ejection of the giant's core from 
the extended envelope, which then rapidly disperses, or by the formation of a common envelope 
binary, which stirs the envelope and leads to its evaporation on a time scale much shorter than the 



giant phase (Livne & Tuchman |1988| ; Davies & Benz |1991| ; Rasio & Shapiro |199Q| , [19911) . These 
studies concentrated on low mass ~ IMq giants and on low velocities that are typical of globular 
clusters, so it is not quite clear how the results scale to higher giant masses and velocities which 



are typical of the GC. More recently, Davies et al. ( 1998 ) proposed that 3-body collisions between 
binaries and giants may provide a more efficient way of destroying the envelope. They specifically 
tailored their dynamical simulations to the velocity field and late-type giant masses and envelope 
sizes that are typical of the GC, and found that envelope-destroying binary-giant collisions require 
Xe of order unity. They also concluded that the role of such collisions is negligible in the fiat core of 
the fiattened isothermal density model, even if all the dwarf stars are in binaries. These works did 
not discuss collisional destruction of early-type giant envelopes. However, these are likely to be less 
effective since it appears that a necessary requirement is that the target star should be composed of 
two distinct components, a compact core and an extended envelope with a shallow density profile. 



^ The K luminosity of a star of radius R increases with the effective temperature roughly as TcH, so that the 
K luminosity scales as B?Tc«, while the bolometric luminosity scales as R^T^g. A giant with R — 100 i?© and 
Teff = 3 10^ K that is reduced by a collision to ~ 1 Rq will maintain its bolometric luminosity by raising its effective 
temperature to Tcff — 3 10'' K, but its K luminosity will fall by a factor of 10^ (7.5™). 
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In the following analysis we will assume that 2-body collisions result in the permanent removal 
of the envelope of a late-type giant if Xg < 0.25. We estimate the efficiency of 3-body collisions 
by assuming that half of the stars are in binaries. If envelope-destroying binary-giant collisions 
typically require = 1, then the effective impact parameter per dwarf star is = 0.5 (1 binary 
per 4 stars, with the collisional cross-section scaling as x^) . The effective impact parameter per 
dwarf star must lie somewhere between the 2-body and 3-body values. Therefore, by modeling both 
types of collisions as 2-body collisions, we are making a conservative estimate of their effectiveness. 

Consider first a star on a circular orbit of radius r around the BH that has total lifetime T and 
becomes a giant with a compact core / extended envelope structure after a time tg. The Poisson 
probability for a giant to survive envelope-destroying collisions with dwarf stars up to time i > is 



S{t) = exp qiT)dTj , (14) 

where q is the local collision rate. The collision rate is the sum of a geometrical cross-section term 
and a gravitational focusing term (see e.g. Binney Sz Tremaine |1987D , 

where Up is the number density of the dwarf projectiles, a is the 1-dimensional velocity dispersion, 
and Mj, Mp, Rt and Rp are the masses and radii of the giant targets and dwarf star projectiles, 
respectively. Both terms have non-negligible contributions over the wide range of stellar masses. 



radii and velocities in the GC (Eq. |B^). Equation |15| assumes a Maxwellian DF with an isotropic 
velocity dispersion. These conditions appear to hold in the inner GC. The observed radial velocity 
distribution is well fitted by a Gaussian distribution (Genzel et al. 1996| ) and is thus consistent 



with theoretical predictions of the velocity distribution very near a BH (Quinlan, Hernquist & 
Sigurdsson 1995| ). The measured velocities in the plane of the sky and along the line-of-sight are 



consistent with an isotropic velocity dispersion (Chez et al 1998) 



The collision rate q depends on the position r through Hp and a. The observed isotropic 
velocity dispersion in the inner GC implies that stars observed at r have a wide range of orbital 
parameters and spend a large fraction of their orbital period at distances much larger or much 
smaller than r. It is therefore necessary to take into account the fact that the collision rates are 
not constant along the stellar orbits. As we show in appendix |^, when the density and velocity 
dispersion are power-laws in r, the orbital averaging modifies Eq. ^ only by introducing two 
constant correction factors for each of the geometrical cross-section term and the gravitational 
focusing term. These factors depend on the power-law indices and on the ratio ^ = v'^/a'^, where 
Wg = GM,/r is the circular velocity near the BH. 

In our analysis we assume that stars that are not bound to the BH are unaffected by collisions, 
since the transit time through the inner 0.1 pc is only ~ 1000 yr, which is much shorter than the 



- 16 - 



typical mean time between collisions (> 10^ yr). Therefore, the orbit-averaged effective collision 
rate depends very strongly on the value of ^. The larger this ratio is, the smaller the fraction of 
unbound stars, the more tightly bound the stars whose orbits pass through r and the longer they 
spend in the high density regions of the inner GC. The ratio ^ is related to the stellar density by 
the Jeans equation, which for a steady-state, non-rotating system with an isotropic potential and 
isotropic velocity dispersion is 

^ v1 dhin dlna"^ , , 

^ = ^ = -T^ In • 16 

cr^ a In r a in r 

Near the BH v'^ ~ GM,/r, so the solution of the Jeans equation for n ~ r^" is cr^ ~ 1/r and 
therefore £, ~ 1 + a. 



Genzel et al. (|1996 ) assumed a flattened power-law density model (§ |2.3.3|) , which is flat (a ~ 0) 



near the BH, and used the Jeans equation to estimate ^ from the observed projected velocity 
dispersion ap. They assumed a deprojected velocity dispersion of the form cr^ = cj^ -|- a'^{r/rc)^'^ 
with (Too a fixed constant, and obtained a best fit value of 7 = 1.2 and consequently a ^ = v'l/a'^ 
ratio that falls from ~ 2 at r ~ Tc down to 1.2 at r = 0. However, their velocity dispersion 
model is inconsistent with a density cusp. For example, the Jeans equation associates a n oc 
cusp with a smaller deprojected velocity dispersion that varies as o"^ = GM,/3r. This can be 
understood qualitatively by noting that in a flat core the stars have a wider distribution along the 
line-of-sight and are on average farther away from the BH than they are in a cusp. Therefore, in a 
flat core a large value of Up implies that the stars move at a large fraction of the circular velocity, 
whereas in a cusp this merely reflects the fact that a larger fraction of the stars are deep in the 
potential well. A reduction of the velocity dispersion significantly changes the orbital structure of 
the population. The fraction of stars on orbits more bound than circular {v^ < ^c); l^ss bound 
than circular {v^ <v^< 2v1) and unbound (2u^ < f^) are respectively 0.20 : 0.23 : 0.57 for ^ = 1, 
0.43 : 0.31 : 0.26 for ^ = 2, and 0.61 : 0.28 : 0.11 for ^ = 3. That is, the steeper the cusp, the more 
the stars are bound to the black hole. This implies that the stellar population in a flat core is 
significantly more mixed than it would be in a cusp, and that any local spectral variations in the 
stellar population will tend to be diluted. 

The probability that a star observed at a random time t during its lifetime will be brighter 
than some magnitude Kq is given by 



Po = ^yJj mo{t)dt + ^ mo{t)S{t)dtj , (17) 

where mo{t) is a function that equals 1 when the star is brighter than Kq, and equals otherwise. 
The probability Pq depends on the star's evolutionary track through the fraction of time the star 
spends in the giant phase and through the evolution of the stellar radius and the -fC-luminosity. 
The probability Pq has to be averaged over the stellar population to obtain the total probability 
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of observing any star brighter than Kq. As was discussed in the stellar population in the 
GC that contributes to the KLF is well modeled by a continuous star-forming history with an 
approximately Salpeter initial mass function with a minimal mass of Mi = O.SMq and a maximal 
mass of M2 = I2OM0, so that the present-day mass function (PMF) can be approximated as 

-ito)oc^i ,.-2.35. ' (18) 



where ^ is the star-formation rate, which we take to be constant, Mq is the zero-age main sequence 
stellar mass, and to is the age of the system, which we take to be 10 Gyr. 

We calculate the PMF-averaged Pq using the Geneva stellar evolution tracks (Schaerer et al. 



19931) for twice-solar metallicity stars and assume blackbody stellar spectra for tracking the K 
luminosity with time. We model the stellar DF by a two power-law function (Eq. |6|) with = 10" 
and /3ft = (3 — a) X IO^Mqpc"^, where the pre-factor normalizes all models to have an equal mass 
of ~ 8 X 10^ Mq within r^. For the dwarf projectiles (burning stars and stellar remnants), we 
assume np{rh) = Pb/ ("i) with (m) = 0.6Mq, while for the target stars (burning stars only), we 
assume nt{rh) = pb/A (see Alexander & Sternberg 1999| ). The number of stars per p-annulus that 
are brighter than Kq is then given by 



/+00 f-M'z 
/ nt{p, z) [Po(l - fub) + fub] Ttr^^^o^^ ' (19) 

where fub is the fraction of unbound stars in the population and where Pq is calculated with the 
effective orbit-averaged collision rate as described in appendix 

We begin by considering a cusp (^ = 3 and fuh = 0.11). Figure ^ overlays the distribution 
of the observed stars on a binned gray-scale plot of N{p;Kq). The middle contour line shows the 
model prediction for the values of the minimal magnitude above which less than one star per p-bin 
is expected to be observed, with and without collisions. The n oc cusp model with collisions 
reproduces the observed trend to a remarkable degree, given that no attempt was made to fit the 
data. Apart for the Hel/Hl stars, which clearly constitute a separate sub-population, and the 
other young stars (early- type and dust-embedded stars), there are in total 24 stars below the '1 
star per p-bin' contour. This is statistically consistent with the expected value of 18 for the 18 
p-bins between 0" and 4.5", especially given the possibility that a few of the 24 stars could be 
unidentified early type stars. The minimal magnitude contour for the no-collisions case clearly fails 
to reproduce the observations, as it remains flat down to the smallest projected radii. Our model 
for the collisional depletion of giants predicts that the collisions become effective only at p ^ 2", 
which corresponds to > 4 x 10''pc~^ and a > 300 km s^^. 

Up to this point we considered only the collisionless orbits of the stars. However, two-body 
encounters will also affect the orbits. The typical mean time between envelope-destroying collisions 
in the inner 2" is ~ 10^ yr per giant, whereas the crossing time of the collisionally-depleted region 
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is only ~ 500 yr. This raises the question whether the central region can avoid being replenished 
by collisionally-deflected giants from the dense stellar environment just outside it. A detailed 
calculation of the replenishment rate is beyond the scope of this work, but we note that there are 
several reasons that suggest that this process is likely to be inefficient. Slow diffusion of giants 
into the inner region may make them as vulnerable to collisional destruction as the giants that 
were originally bound to that region. Giants that will be deflected by a single close elastic collision 
into a low angular momentum orbit will spend only a small fraction of their lifetime in the central 
2". Moreover, the velocity dispersion in the cusp is high, so that large deflection angles require 
very close encounters. This defines a collisional radius, TcoII (Frank & Rees 1976| ), inside which 
large defiections require physical (inelastic) collisions. For a Rt = WO Rq, Mt = 3Mq giant, the 
collisional radius is 



rcoii = ^^t ~2pc~50". (20) 

Therefore, a close encounter within Xgrcou, which is close enough to defiect the giant into a low 
angular momentum orbit, will also destroy its envelope. Replenishment by giants from beyond 
iCe'^coU will be less effective because of the lower stellar density there and because of the small 
angular size of the central 2 , as seen from v > x^r^QH. 

Next, we consider the case of a fiat core (a = 0, /? = 2, = 10" with ^ = 1 and fub = 0.57). 
Surprisingly, the fiat core model shows that the apparent depletion of the high luminosity stars 
can also be a projection effect. Figure |9| shows that collisions play no role in such a low density 
models, in agreement with the conclusion of Davies et al. ( |1998 ) but contrary to the suggestion of 



Genzel et al. (|1996| ). However, since the number of stars per p-annulus falls as p in a fiat core, the 
probability of observing one of the rare luminous star also falls with p. This is to be contrasted 
with the cusp, where the number of stars per annulus is constant. 

The origin of the depletion in a cusp is collisional, while in a fiat core it is due to 
a projection effect^. Stellar cusps shallower than will display both types of depletion. 
Close inspection of the predicted depletion in the cusp (Fig. |9|) shows that it somewhat 
underestimates the depletion at the smallest projected distances, and that a better fit can be 
obtained by a shallower cusp, for example r~^/^. 

We conclude that the depletion of the bright stars in the inner 2" can be explained by a 
combination of collisional destruction of giant envelopes, which requires the high densities of a 
stellar cusp, and by projection effects. The latter is due to the fact that density distributions that 
are shallower than a cusp have progressively less stars enclosed in smaller projected radii, and 
so there is a smaller chance of observing the rarer, more luminous stars there. Because of the 
the degeneracy between the collisional depletion and the projection effects, the distribution of the 



A r " cusp with a > 2 will have collisional depletion offset by apparent enhancement due to projection. 
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brightest stars alone cannot decide the case between a flat core and a cusp. 



6. Discussion 

6.1. The star counts 

The analysis of the star counts is not straightforward and must overcome several difficulties. 
First, the observed star-counts have yet to be published in a way suitable for detailed analysis by 
taking into account completeness corrections and field edge corrections. Second, the GC contains 
a population of luminous young stars, which are not dynamically relaxed, and therefore introduce 
noise in the counts. Third, a gradient in the extinction field may bias the counts. Although none of 
these concerns can be completely eliminated, we have attempted to minimize their possible effects 
in our analysis. 

An important aspect of our effort to control the effects of incompleteness was the use of three 
different data sets, with different properties and biases. We analyzed them in a uniform way by 
methods that avoid binning, as is appropriate when dealing with sparse data. It is reassuring 
that all three data sets indicate that a cusp matches the counts better than a flat core. The case 
for a cusp is further strengthened by the fact that both incompleteness due to source confusion 
and a central concentration of dust will tend to suppress the observed cusp. Another way of 
circumventing the completeness problem is to consider the information in the distribution of the 
brightest stars. We discuss this further below. 

We attempted to limit the contamination by the young unrelaxed population by excluding all 
the identified young stars from the data sets. Although we do not know how many of the fainter 
stars are associated with this unrelaxed population, the star-burst model of Krabbe et al. ( 1995[) 



predicts that their fraction should become negligible as the counts go deeper. It is possible that 
this trend is already apparent in the differences between the OSIRIS data set and the two other 
deeper data sets (Figure]^. The ML power-law index does not change significantly when the 
young stars are included in the sample (Figure]^. 

To estimate the possible effects of extinction gradients, we constructed a few simple 
3-dimensional dust distribution models and tested them against the observed distribution of 
stellar extinction in the projected distance-reddening plane. We find that on the ~ 10" scale, the 
observations are consistent with an external dust screen, which is unlikely to significantly bias the 
derived power-law index of the stellar distribution. Our constraints on variations on ~ 1" scales 
are weaker, but we note that the detection of moderately reddened stars very deep in the potential 
well is consistent with little reddening in the innermost GC, and that no significant color gradient 
was detected in the inner 1". 



These precautions and checks are no substitute for better data. Clearly, deep, completeness- 
corrected data sets are desirable. A comparison of the likelihood curves of the Keck and 
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3D+SHARP data set shows that even a two-fold increase in the sample size from ~ 50 to 100 
can lead to a significant improvement in the ability to discriminate between stellar distribution 
models. Our KLF model suggests a two-fold increase in the number of stars observed for every V 
improvement in the photometric sensitivity. 



6.2. The depletion of bright stars 

In addition to the analysis of the total star counts, we examined the distribution of the counts 
as function of magnitude and noted that the highest resolution data set shows a gradual depletion 
of the brightest stars in the inner 2" (a reflection of the fact that there are no bright stars in 
the Sgr A* cluster). Reproducing this effect by reddening would require a ~ 2.5™/2" extinction 
gradient, which is not indicated by our analysis (§^. We showed that envelope-destroying collisions 
between dwarfs and giants can reproduce the observed trend to a remarkable degree, once the 
young unrelaxed stars are excluded from the data set. Our detailed model for the effects of 
collisions on the luminosity function included six components: realistically lowered estimates for the 
cross-section of envelope-destroying collisions, based on results from hydrodynamical simulations 
of dwarf-giant collisions; an orbit-averaged estimate of the collision rate; a velocity dispersion 
model that is self-consistent with the assumed density distribution; a stellar mass function that 
reproduces the mean observed K luminosity function and mass-to-light ratio; stellar tracks that 
model the evolution of the size and luminosity of the stars as function of their initial mass; and a 
careful statistical treatment of the survival probability. 

Our dynamical analysis shows that a cusp is necessary for collisional depletion to play a role 
because the process becomes effective only at a density of order > 5 x IO^Mqpc^^ with a velocity 
dispersion of order > 300kms~^ , and because a cusp is much more tightly bound to the BH 
than a flat core. The unbound stars spend too little time in the high density central regions to be 
affected by the collisions, and so set an upper limit on the fraction of giants in the population that 
can be effectively destroyed. 

Stellar collisions are only one way to explain the absence of bright stars without having to 
invoke intrinsic variations in the stellar population. We also find that a flat core can qualitatively 
mimic such a depletion by a projection effect, even when collisions play no role. This explanation 
for the absence of bright stars in a flat core may actually be more likely than to assume intrinsic 
variations in the stellar population. Any localized, distinct stellar population in a flat core will be 
diluted by the large fraction of loosely bound stars that sample the population from over the entire 
GC. 
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6.3. Implications of a stellar cusp in the GC 



We now discuss briefly some of the implications of the existence of a steep stellar cusp in 
the GC. The tidal disruption rate of stars by the BH is enhanced in a steep cusp. Such events 
may be observed as short and luminous flares (Rees 198^ ). Bahcall & Wolf (|1977| ) estimated the 



diff'usion rate into the 'loss-cone' (low angular-momentum orbits) as function of the BH mass 
and of the line-of sight velocity dispersion and the stellar density outside the cusp. Adopting a 
broken power-law density distribution (Eq. with = 10", (3 = 2 and ni, = 10^ pc~^ and taking 
Up ~ 85kms~^ at ~ 2.5 pc from Kent ( |1992| ) or dp ~ 60kms~^ at ~ 4 pc from Genzel et al. 
(1996), we estimate a loss-cone diffusion rate of few x 10~^yr~^. This is consistent with the rates 



derived by Magorrian & Tremaine ( 1999| ) for faint galaxies {L < 10^^ Lq) with low-mass BHs and 



steep density cusps, but is much larger than the 10~^ — 10~^yr~^ rates that they derived for L* 
galaxies. Magorrian & Tremaine suggest that tidal flares are most likely to be detected in faint 
galaxies. Our results raise the possibility that at least some luminous galaxies may have tidal flare 
rates as high as those of faint galaxies. 

Inelastic dwarf-dwarf collisions are expected to modify the stellar distribution in the inner, 
highest density regions of the cusp by gradually destroying the tightly bound stars. This leads to 
the formation of a n oc r~^/^ cusp of unbound stars. The relevant time-scale for this process is the 
time it takes to destroy a star by successive collisions. Murphy, Cohn and Durisen ( [1991| ) show that 
the collisional destruction timescale is ^ lOicoU, where tcoU is the collision timescale per star. The 
size of the unbound cusp that they derive in two of their numerical models that roughly bracket the 
conditions in the GC (models 3B and 4B), is 0.002 < < 0.04 pc at 15 Gyr. This qualitatively 
agrees with our independent estimate of icoU ^ 5 x 10* yr at r < 0.02 for dwarf-dwarf collisions in 
a a = 1.5 cusp. We therefore expect the cusp to flatten somewhere within p ^ 0.5" . We estimate 
that dwarf-dwarf inelastic collisions occur in that volume at a rate of 10~^ — 10~^yr~^. 

A steep stellar cusp will also increase the rate of microlensing of stars by the BH (Alexander 
& Sternberg 1999|) , and especially that of short duration microlensing events of stars very near the 



BH. Such microlensing light-curves can be used, in principle, to probe the mass distribution of the 
central dark mass and test whether it is a point mass. These issues will be discussed in more detail 
elsewhere. 

Finally, we note that the integrated stellar mass close to the black hole is very small in spite 
of the very high density of the cusp. For example, the stellar mass within 2" in a cusp is 
~ 8 X 10^ Mq as compared to only ~ 8 x 10^ Mq in the flattened isothermal distribution (when 
both are normalized to have the same mass within Vc = 10"). This difference amounts to only a 
~ 3% over-estimate in the BH mass due to the fact that some of the stellar mass is erroneously 
attributed to the BH. 
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7. Summary 

Recent deep high-resolution observations of the inner GC and large scale multi-color 
observations have supplied a wealth of information on the spatial distribution of stars near the 
super-massive BH in the GC. This offers the opportunity to map the stellar distribution around 
the BH directly from the star counts, rather than from the integrated light. Dynamical scenarios 
predict that a stellar cusp should form around the BH, and that the shape of the cusp could 
provide clues on the formation history of the BH and the stellar system around it. It is therefore 
puzzling that recent empirical GC mass models indicate a flat core rather than a cusp, while at 
the same time, deep observations reveal an over-density of faint stars and an absence of bright 
stars in the inner 1" (the Sgr A* cluster). Taken at face value, this picture raises several questions: 
What does the absence of a cusp imply about the dynamical state of the GC? Are the faint stars a 
distinct population? If so, how have they formed in the inner l" or how have they migrated there? 
Our analysis of three independent sets of star counts points to a simpler interpretation of these 
findings. Our main results are as follows. 



1. The central over-density is likely the tip of a stellar cusp that rises smoothly in the inner 
~ 10". A stellar cusp is consistent with the star-counts and is preferred at the ~ 2a level 
over a flat core distribution even by the least restrictive data set. A flat core is completely 
ruled out by the largest data set. The maximum likelihood range for the power-law index 
of the cusp is a ~ 1.5 to 1.75. An alternative way of stating this result is that our analysis 
points to a small core radius, Vc < 2.5" . 

2. The gradual depletion of the luminous giants in the inner 2" can be explained by two effects: 
envelope-destroying dwarf-giant collisions in a steep, dense and tightly bound cusp, or a 
projection effect in a flatter distribution, where the number of stars per projected annulus 
falls with decreasing projected radius. The observed depletion in the innermost regions is 
better described by a a < 2 cusp, where both effects contribute. Either way, it is unnecessary 
to assume that the stellar population in the innermost GC is intrinsically different (fainter) 
from the general population. 

3. The joint constraints of the star counts and the depletion of the brightest stars suggest that 
a power-law cusp with a in the range 3/2 to 7/4 provides a good overall description of the 



data. Such a cusp is consistent with the Bahcall-Wolf solution (Bahcall & Wolf |1976| , |l977l ; 



Murphy, Cohn & Durisen |1991| ; see § p.lD for the distribution of intermediate mass stars in 
a multi-mass stellar system that has undergone two-body relaxation around a black hole 
(Eq. |2|). Dynamical estimates (Eq. ^ suggest that the inner GC is relaxed. If future deeper 
star counts support these conclusions, this would be the first detection of such a system. 



I am grateful to A. Eckart, R. Genzel and A. Chez for supplying additional information about 
their data and to J. Bahcall, D. Eisenstein, J. Goodman and S. Tremaine for helpful comments. 
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A. The stellar "ring" at 2" 

Is the striking over-concentration of stars in a thin ring at p ~ 2", which is seen in the 
cumulative distribution of the Keck data set (Figs. |l] and inconsistent with a power-law 
distribution (9 late-type stars in a ring of width 0.06" out of a total of 54 late- type star or 11 stars 
of all spectral types out of a total of 64 stars)? This is probably not an instrumental artifact, since 
it is seen also in the full OSIRIS data set. The 11 stars do not seem to have any other peculiar 
characteristics apart from their spatial alignment. Only two (16CC and 29N) are identified as 



Hei/Hi stars (Eckart & Genzel |1997|) , their mean magnitude of K = 12.6™ ±2.9™ is consistent with 
the total average of K = 13.4™ it 2.8™, and their projected velocity vectors appear to be random. 
In order to address this question, we performed Monte-Carlo simulations where we randomly drew 
64 stars from a n oc r'" distribution and recorded the number of times a group of k stars or more 
was concentrated within a ring of Ap = 0.06" or less, for various values of k and a (a ring of k 
stars is counted twice if it is a subset of a ring of k + 1 stars). The distribution of the ring radii 
was also recorded. The results, listed in table |2|, show that a ring of 11 stars in 0.06" is quite rare 
in a a = 2 model, where there are on average only 2.3 10"^ such rings per random realization, with 
radii 1.3" ± 0.7", but less rare in the flatter a = 3/2 model, where the average number is 2.3 x 10~^ 
with radii 2.1" it 0.3". The probabilities increase by more than an order of magnitude if only the 9 
late-type stars are considered members in the ring. 

The Monte-Carlo averaged ring position reflects a trade-off between the ring area, which 
increases with p, and the density distribution, which decreases with p. Therefore, both the 
probability for observing a ring, and the ring's radius, increase in flatter distributions. The fact 
that an 11-star ring of radius 2" is observed, can be turned around and used as a Maximum 
Likelihood argument in favor of a flatter inner distribution. Future deeper number counts will 
show whether this feature is more than a statistical fluctuation. 



B. The effective orbit-averaged collision rate 

We estimate the corrections to the local collision rate (Eq. |l5|) due to the fact that the stars 
are not on circular orbits by making two simplifying assumptions. First, we assume that the stars 
in the inner GC move on Keplerian orbits around the BH, so that the orbits can be parameterized 
by their semi-major axis a and their eccentricity e. Second, we assume that unbound stars are 
unaffected by collisions (see When the orbital period P of a star around the BH is shorter 
than the giant phase, the orbit-averaged collision rate is 

2 ra{l+e) 

q{r,a,e) = q{s)— , (Bl) 

^ Ja{l~e) S 
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where s is the radial coordinate of the orbit and a(l ± e) are the apoastron and periastron radii, 
respectively. The radial velocity varies with s as 



2 e + 



GM\ 
s J s"^ 



(B2) 



where e = —GM,/2a is the specific energy and h = \J GM,a(l — e^) is the specific angular 
momentum of the star. The orbital parameters of a star at r with velocity v can be expressed in 
terms of v = v/cj and the angle between the radius and the velocity vectors, t/, 

(B3) 



1 + sin T] 



2f_ 



where ^ = vl/a"^. For n oc r~" and cr^ oc r"''', the geometrical cross-section term 
and the gravitational focusing term in the local collision rate can be expressed as 
q{r) = Qg + Qf = Ar~^'^~^'^^'^^ + Br~^°'~'^/'^\ where the constants A and B are defined by 
Eq. In this case, the corresponding orbit-averaged terms are 



vr 



X2 ^a±7/2-2^^ 
XX \/{x- Xl)(x2 - x) 



(B4) 



= Cgjiv, i])qgj{r) , (B5) 



where xi^2 = 1/(1 i e) and where the (+) case corresponds to qg and the (— ) case to qf. It then 
follows that {q)g j, the averages of qgj over the Maxwellian DF of the bound stars, are related to 
the local collision rates qgj by proportionality factors that are functions of the constants a, 7 and 
^ but not of r. 



^^^9J ^ I ^9 J = (2vr)-'/' ^ Cgj{v,ii)v^ exp {-v^ /2) sin r? dvdr^ . (B6) 

In practice, the expectation value (C) diverges for certain choices of a and 7 due to low 
angular-momentum orbits that sample regions of diverging density and velocity dispersion. For 
this reason, and because the survival probability (Eq. [ij) is a highly non-linear function of the 
collision rate, a more relevant quantity is the median correction factor C, defined as the ratio 
between the median of q over the population of bound stars and the local rate q. Table ^ lists 
Cf and Cg for different density distributions that may be relevant in the inner GC, as well as the 
fraction of unbound stars, fuh-, that corresponds to each value of ^. The median correction factors 
were estimated by 10® Monte Carlo draws of random velocity vectors from the DF. Using these 
factors, the effective collision rate for the bound stars is estimated by 



Qeff 



ACgT- +BCfr- . (B7) 
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The relative contributions of the geometrical cross-section term and the gravitational focusing term 
to the effective collision rate is (cf. Eq. ^) 



qg\ _Cg2 [ M, \ (xeRt + Rp 



'ifJeS Cfi\Mt + Mp 



(B8) 



This ratio can be both much larger than 1 or smaller than 1 over the wide range of stellar masses, 
radii and distances in the GC, so that both terms should be taken into account when calculating 
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Fig. 1. — 90 stars in the inner ~ 6" x 6" observed by the Keck interferometer (Ghez et al. 1998 ). 
Stars marked by filled circles constitute the "ring" feature (Appendix §^. Only the 64 stars in the 
central 2.5" (of which 10 are early-type stars) are used in the analysis. 
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Fig. 2.— The inner 40" x 40" of the OSIRIS field (Blum et al. [T996| ). Stars marked by filled circles 
are highly reddened stars (Ak > 4.5™). Only the 50 late-type stars in the inner 13" are used in the 
analysis. 
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Fig. 3.— The combined 3D (Genzel et al. [1996| ) and SHARP (Eckart & Genzel [19971 ) counts in 
the inner 20". Stars marked by filled circles are identified Hei/Hi and early- type stars. Only the 
92 late-type 3D stars and the 15 spectrally unidentified SHARP stars (excluding the faint SgiA* 
cluster stars) in the inner 13" are used in the analysis. 
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Fig. 4. — The cumulative DF of the late-type stars for the different data sets, compared to various 
stellar density distribution models. Top left: the inner 2.5" of the Keck field (Ghez et al. 1998 ). 
Top right: the inner 13" of the OSIRIS field (Blum et al. |1996D . Bottom: the inner 13" of the 3D 
field (Genzel et al. |1996| ) together with the inner 10" of the SHARP data excluding the faint Sgr A* 
cluster stars (Eckart & Genzel 1997| ). 
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Fig. 5. — The ML estimate of the stellar distribution for the late-type stars in the three data 
sets. The likelihood curves for the full Keck data set (including early-type stars) is also presented 
for comparison. Left: The likelihood curves for the power-law index a of the broken power-law 
distribution, for fixed values of /3 = 2, = 10". Right: The likelihood curves for the core radius Tc 
of the flattened isothermal distribution. 
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Fig. 6. — The extinction of individual stars in the inner 25" of the GC deduced from the near 
infrared colors (Blum et al. 1996| ). The average error on the extinction estimate is ~ 0.1™ (not 
including systematic error due to the assumed intrinsic colors). The linear regression of Ak on p in 
the inner 10" and 20" is also shown. IRS 21 is a highly polarized extended object (Ott, Eckart & 
Genzel 1199^ ). 
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Fig. 7. — A simple trial model for the extinction field in the GC. A star at projected distance p 
and point z along the line-of-sight, is embedded in a spherical dust cloud of core radius R that 
surrounds the black hole. The total extinction is the sum of the extinction due to dust inside the 
cloud, A2{p,z), and to dust in a foreground screen, Ai{p). 
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Fig. 8. — The relative probability density (normalized to the peak) for the detection of stars at 
projected distance p and extinction Ak- The dark regions correspond to high detection probability. 
Top: a n oc stellar distribution. Bottom: a flattened isothermal model (Eq. ^). The assumed 
dust model is a spherical d = 1 / (l + 2,{r/Rf) distribution with R = 30", D = 0.15, and Ai = 2.5°", 
the luminosity function is df /dK oc lO''^ with h = 0.35, and the detection threshold is Kq = 13™. 



Fig. 9. — The mean number of stars with K < Kq per 0.25"-wide bin, as expected with colhsional 
destruction of giants (gray-scale and full line contours) and without collisional destruction (dashed 
line contours). The contour lines mark, from top to bottom, the model predictions for the minimal 
magnitude below which less than 1.5, 1 and 0.5 stars per 0.25" bin are expected to be observed, on 
average. Top panel: a stellar distribution. Middle panel: a flat core in the inner 10" and a 
distribution outside. Bottom panel: a r~^/^ distribution in the inner 10" and a distribution 
outside. Overlayed on the gray-scale are the stars observed by Ghez et al. ( 1998|) , marked by 



different symbols according to their spectral type (Eckart & Genzel |1997| ; Genzel et al. [19971 ): 
Hei/Hi stars (H), late-type stars (L), early-type stars (E) , dust embedded stars (D), and type 
unknown (circle). 
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Table 1: The Kolmogorov-Smirnov acceptance probabilities for various stellar density models. The 
broken power-law models have ri, = 10" and the flattened isothermal model Vc = 10". The highest 
probabilities for each data set are emphasized in italics. 



data set 


Keck 


Keck 


OSIRIS 


3D+SHARP 


Subset 


All 


old stars 


old stars 


old stars w/o Sgr A* 


No. of stars 


64 


54 


50 


107 


Flat, isothermal 


0.05 


0.12 


0.41 


0.20 


r-2 


0.02 


0.02 


0.10 


0.04 




0.11 


0.12 


0.78 


0.63 




0.38 


0.31 


0.33 


0.50 


fp 7/4 _|_ J, 2 


0.09 


0.10 


0.51 


0.34 


rp 3/2 _|_ J, 2 


0.27 


0.27 


0.91 


0.56 



Table 2: The frequency of ring-like structures in 10 Monte-Carlo realizations drawn from a n oc r~" 
density cusp. In each realization, 64 stars are randomly distributed within 2.5", and the number of 
rings of width < 0.06" is recorded as function of the number of stars in the ring. 



Stars per ring Rings per realization Mean ring position 



a = 3/2 


> 7 


1.8 10"^ 


2.0" ± 0.4" 




> 8 


4.010~2 


2.0" ± 0.4" 




> 9 


8.110-3 


2.1" ±0.3" 




> 10 


1.510-3 


2.1" ±0.3" 




> 11 


2.310-^ 


2.1" ±0.3" 


a = 2 










> 7 


5.410-2 


1.3" ± 0.7" 




> 8 


9.210-3 


1.3" ± 0.7" 




> 9 


1.410-3 


1.3" ± 0.7" 




> 10 


1.910-^ 


1.3" ± 0.7" 




> 11 


2.310-^ 


1.3" ± 0.7" 



Table 3: The orbital correction factors to the local collisional rate (Eq. p7| ) and the fraction of 
unbound stars in the population. 



a 


7 




Cf 


^ = a ± 7 


fub 


5/2 


1 


6.2 


2.3 


7/2 


0.07 


2 


1 


3.5 


1.7 


3 


0.11 


7/4 


1 


2.8 


1.5 


11/4 


0.14 


3/2 


1 


2.2 


1.4 


5/2 


0.17 


1 


1 


1.6 


1.2 


2 


0.26 


1/2 


1 


1.1 


1.1 


3/2 


0.39 





1 


0.9 


1.0 


1 


0.57 



